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ABSTRACT 

We report the detection of HCO + (J=1^0) emission towards the Cloverleaf quasar (z = 2.56) 
through observations with the Very Large Array. This is the first detection of ionized molecular gas 
emission at high redshift (z > 2). HCO + emission is a star formation indicator similar to HCN, 
tracing dense molecular hydrogen gas (n(H2) ~ 10 5 cm -3 ) within star-forming molecular clouds. We 
derive a lensing-corrected HCO + line luminosity of L' aco+ = 3.5 x 10 9 K kms _1 pc 2 . Combining our 
new results with CO and HCN measurements from the literature, we find a HCO + /CO luminosity 
ratio of 0.08 and a HCO + /HCN luminosity ratio of 0.8. These ratios fall within the scatter of the 
same relationships found for low-z star-forming galaxies. However, a HCO + /HCN luminosity ratio 
close to unity would not be expected for the Cloverleaf if the recently suggested relation between 
this ratio and the far-infrared luminosity were to hold. We conclude that a ratio between HCO + 
and HCN luminosity close to 1 is likely due to the fact that the emission from both lines is optically 
thick and thermalized and emerges from dense regions of similar volumes. The CO, HCN and HCO + 
luminosities suggest that the Cloverleaf is a composite AGN-starburst system, in agreement with 
the previous finding that about 20% of the total infrared luminosity in this system results from dust 
heated by star formation rather than heating by the AGN. We conclude that HCO + is potentially a 
good tracer for dense molecular gas at high redshift. 

Subject headings: galaxies: active, starburst, formation, high redshift — cosmology: observations — 
radio lines: galaxies 



1. INTRODUCTION 

One important goal in studies of galaxy formation is to 
determine star formation characteristics through obser- 
vations of molecular gas in the early universe. Molecular 
gas in high-redshift galaxies is commonly traced by CO 
emission and has been found in >30 galaxies at z > 2 
to date. The observed molecular gas with masses of 
> 10 10 M© provides the requisite material for st ar for- 
mation (see review by Solomon & Vanden Bout 120051 
and references therein). 

CO is a good indicator for the total molecular gas con- 
tent of a system, as it can be excited at relatively low 
densities; its low dipole moment of = 0.1 implies 
a critical density of only nn 2 ~ 10 2 — 10 3 cm~ 3 for the 
lower- J transitions. Hence, emission from low-J CO 
transitions is a relatively poor tracer of the denser gas 
that is more intimately associated with star formation. 
The most common tracers of the dense molecular gas 
phase are HCN and HCO + . Both molecules have much 
higher dipole moments (>g CN = 2.98, ^g co+ = 4.48) 
than CO. The critical density to collisionally thermalize 
their lower-J transitions is therefore much higher than 
for CO, ?ih 2 ~ 10 5 - 10 6 cm~ 3 (e.g. Gao & Solomon 
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I2004all27)04ll Brouillet et al joog) . 

Recent studies of the dense molecular gas phase in 
nearby (z < 0.3) luminous and ultra-luminous in- 
frared galaxies (LIRGs/ULIRGs) have shown that the 
HCN and HCO + luminosities correlate with the star- 
formation rate as traced by the far-infrared (FIR) lumi- 
nosity. These correlations are tighter than the correlation 
betwee n the CO and FIR luminositi es (Solomon et al. 
1992a; Gao & Solomon l2004aL I2004bt Gracia-Carpio et 
al.|2006). Due to the relative faintness of the emission 
lines, HCN has only been detected in four objects at high 
z to date, and all detec tions are in the host gala xies of 
quasars (So lomon et al. 120031 : Vand en Bout et al. 120041: 
Carilli et al. 120051 Wagg et al. l27)05h . 

It has recently been argued that AGN-dominated 
galaxies have higher HCN/HCO+ and HCN/CO lumi- 
nosity ratios than s tarbu rst-dominated g alaxie s (Kohno 
et al. [2001; Kohno 120051: Imanishi et al. 12006ft . In this 
context, it has been suggested that the presence of X-ray 
emission emerging from a dust-enshrouded AGN may 
significantly enhance the chemical ab undan ce of HCN 
rel ative to HCO+ ( Lepp fc Dal garno 1996; Kohno et 
al. 120011 Usero et al. 1200 4). Also, the excitation of the 
HCN molecule may be affected by IR-pu mping through 
a 14 fim vibrational band (Aalto et al. 1995;. Based 
on these considerations and their H CO + s urvey of low- 
z (U)LIRGs, Gracia-Carpio et al. (2006) suggest that 
the HCO+-to-HCN intensity ratio towards FIR-bright 
(Lfir > 10 12 L Q ) objects (like the Cloverleaf) are likely 
low. Measurements of HCO + emission in high redshift 
objects would thus lead to new constraints on their dense 
interstellar medium and, potentially, on the radiation 
field pervading it. 

In this letter, we report the first high-z detection of 
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HCO + (J=1^0) emission, which was observed towards 
the Clover leaf quasar (z — 2.56) with the Very Large 
Array (VLA) 6 . Due to its strong gravitational magnifica- 
tion (magnification factor /iL = 11, Venturini & Solomon 
2003), the Cloverleaf is the brightest CO source at high 
redshift, and it also exhibi ts bright HCN(J=1^0) emis- 
sion (Solomon et al. I2003|) and emission fro m both C i 
fine structure lines (Weiss etal. 120031 l2005|) . A previ- 
ous search for HCO + (J=4^3) emission in this source 
has been unsuccessful, setting an upper limit of 14mJy 
on the peak line flux density (Wilner etal. [1995). We 
use a standard concordance cosmology throughout, with 
H = 73 km s -1 Mpc -1 n M = 0.24, and n A = 0.72 
(Spergel et al. HM |2006|). 

2. OBSERVATIONS 

We observed the HCO + (J=1^0) transition line 
(z/ rest = 89.1885230 GHz) towards H1413+117 (the 
Cloverleaf quasar) using the VLA in D configura- 
tion between 2005 November 26 and 2006 January 13. 
At the target z of 2.55784, the line is redshifted to 
25.068166GHz (11.96mm). The total integration time 
amounts to 25.5 h. Observations were performed in fast- 
switching mode using the nearby source 14160+13204 
(distance to the Cloverleaf: 1.8°) for secondary ampli- 
tude and phase calibration. Observations were carried 
out under excellent weather conditions with 23 anten- 
nas. The phase stability in all runs was excellent (typ- 
ically <15° rms for the longest baselines). For primary 
flux calibration, 3C286 was observed during each run. 

Two 25 MHz wide intermediate frequency bands (IFs) 
with seven 3.125 MHz channels each were observed si- 
multaneously in the so-called 'spectral line' mode cen- 
tered at the HCO + (J=1^0) line frequency, leading to 
an effective bandwidth of 43.75 MHz (corresponding to 
523 kms -1 at 25.1GHz). In addition, two 50 MHz (cor- 
responding to 598 kms -1 at 25.1 GHz) IFs were observed 
in the so-called 'quasi-continuum' mode at 24.7351 GHz 
and 24.7851 GHz to monitor the continuum at 12 mm. 
We chose a setting with both continuum channels below 
the HCO + (J=1^0) line frequency to avoid the signifi- 
cantly worse system temperatures and locking problems 
of the local oscillator (LO) above 25 GHz. The contin- 
uum was observed for one third of the total time to attain 
the same rms as in the combined line channels. 

For data reduction and analysis, the AIVS package 
was used. The two continuum channels were concate- 
nated in the uv/visibility plane. The HCO + (J=1^0) 
line data cube was generated by subtracting a CLEAN 
component model of the continuum emission from the 
visibility data. All data were mapped using the CLEAN 
algorithm and 'natural' weighting without applying a 
further taper; this results in a synthesized beam of 
4.0 "x 3.0" (~28kpc at z = 2.56). The final rms in the 
combined map is 16 /iJy beam -1 for a 34.375 MHz (cor- 
responding to 411 kms -1 ) channel, and 50^iJy beam -1 
for a 6.25 MHz (75 kms -1 ) channel. 

3. RESULTS 

6 The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 



We have detected emission from the HCO + (J=l— >0) 
transition line towards the Cloverleaf quasar (z = 2.56). 
The source appears to be marginally resolved in both the 
continuum and the HCO + line maps. Two-dimensional 
Gaussian fitting yields a continuum peak flux density 
of 343 ±12 //Jy beam -1 . The continuum-subtracted 
HCO + (J=1^0) line map is shown in Fig. 1. The 
cross indicates the geometrical center position of the 
resolved 12 CO(J=7^6) map at a = 14 h 15 m 46 s .233, 
6 = +11°29'43".50 (Alloin etal. Il997]) . The small off- 
set is likely due to the fact that the lens images are not 
equally bright, i.e. the center of intensity is offset from 
the geometrical center position. The deconvolved source 
size from the Gaussian fit is in good agreement with the 
size of the resolved structure seen in 12 CO( J—7~ >6). The 
line is clearly detected at 8a over a range of 34.375 MHz 
(411 km s -1 ). We derive a line peak flux density of 193 ± 
28 /zJy beam -1 . In Fig. 2, four channel maps (6.25 MHz, 
or 75 kms -1 each) of the central 25 MHz (300 kms -1 ) 
of the HCO + (J=1^0) line are shown. At an rms of 
50^Jybeam -1 , the line is detected at 4<r in the cen- 
tral two channels, and the decline of the line intensity 
towards the line wings is clearly visible in the outer 
channels, as expected. We attribute the small offset be- 
tween the peak positions of channels 2 and 3 to observa- 
tional uncertainties rather than to a real velocity gradi- 
ent. We thus derive a HCO + ( J=l— »0) line luminosity of 
-^HCO+ = 3-5 x 10 9 K kms -1 pc 2 (corrected for gravita- 
tional magnification, /il = 11, and the finite source size 
relative to the synthesized beam, which leads to a 30% 
correction based on the source extension seen in CO, see 
Table and its caption for details). 

We summarize our results in Table 1 together with 
the lin e flux es and luminosities for HCN (Solomon 
et al. l2003h and CO (Weiss et al. l200l . Our 
HCO + (J=1^0) peak flux density corresponds to ~80% 
of the HCN(J=1^0) pe ak flux of 0.24 ± 0.04 mJy de- 
rived by Solomon et al. l)2003|) . This corresponds to a 
-^hco+/^hcn luminosity ratio of 0.8, which is consistent 
with unity within the statistical and systematical uncer- 
tainties. They use an extrapolated continuum peak flux 
at 24.9 GHz of S cont (24 GHz) = 0.26 ± 0.03 mJy. The dif- 
ference to our value may be due to problems with their 
extrapolation or calibration errors. However, it is also 
possible that the continuum of the Cloverleaf is variable 
at 25 GHz. 

4. DISCUSSION 

In the following, we discuss relationships between the 
emission observed in HCO + and other molecules and the 
far-IR continuum for a sample of low-z spiral and star- 
burst galaxies (Nguye n- Q-Rie u et al.[l992; Imanishi et al. 
120011 Gao fe Solom on 120041 . lo w-z (U )LIRGs (Gracia- 
Carpio et al.|2(306; Imanishi et al. | 20 06fl. and the Clover- 
leaf (this work; Solomon et al. [2003; Weiss et al. 12003ft 
as shown in Fig. 3. As only an upper limit exists for the 
12 CO( J= 1^0) line emission in the Cloverleaf (Tsuboi et 
al. H299), we here assume that CO is fully thermalized 
up to the 3^2 transition (i.e., -^co(i-o) = ^co(3-2))- 
We do not discuss effects of differential lensing, which 
could distort the intrinsic luminosity ratios, as models 
indicate similar sizes for mole cular and dust emission in 
the Cloverleaf (Solomon et al. l2003|L 
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Figure 3a: L' HCO+ correlates closely with £fir! a lin- 
ear least squares fit (excluding the Cloverleaf) yields 
log(L Fm ) = (1.11 ± .06) i x log(L'„ rn+ ) + 2.2. For 
HCN, Gao & Solomon ll2004aD find log(L Fm ) = 0.97 x 
log(L^ CN ) + 3.1 based on a larger sample of local star- 
burst and spiral galaxies. Both slopes are consistent with 
unity. 

Figure 3b: £hco+ a ^ so correlates closely with i^ CN ; a 
linear least squares fit (again excluding the Cloverleaf) 
yields log(L^ CN ) = (0.94±0.06)xlog(L^ co+ )+0.6. Pan- 
els a and b of Fig. 3 thus exemplify that HCO 4 " traces 
dense molecular gas as well as HCN, and its close corre- 
lation with the far-IR continuum emission suggests that 
HCO + may also be used as a star formation indicator. 
It is remarkable how well the Cloverleaf agrees with the 
correlations found for local galaxies ranging over more 
than three orders of magnitude in far-IR luminosi ty. 

Figure 3c: Base d on data from K ohno et al. l)200lT) 
and Kohno ((2(305), Imanishi et al. l)2006() argue that 
AGN-dominated galaxies have higher HCN/HCO+ and 
HCN/CO ratios than starburst-dominated galaxies. In 
this diagram, AGN-dominated galaxies thus fall on the 
upper right side, and starburst-dominated galaxies fall 
on the lower left side. For the local starburst/spiral sam- 
ple, we have t aken t he L' nCN / L'^^ rat i° s horn Nguyen- 
Q-Rieu et al . (I1992T) and the L' HCN /L' co ratios from Gao 
& Solomon j2004a) . The Cloverleaf clearly falls on the 
lower left side of the diagram, putting it in the region 
of 'starburst-dominated' galaxies. Indeed, by decompo- 
sition of the dust spectrum into a war m (115 K) and a 
cooler (50 K) component, Weiss et al. !'2()(Kj) find that 
about 60% of the dust emission emerges from the cooler 
component, which may well be dominated by heating 
from star formation. Howev er, based on an Arp220 tem- 
plate, Solomon et al. |2Q03f) have shown that only about 
20% of the total far-IR luminosity is powered by the 
starburst. It is therefore unclear whether or not this di- 
agram should indeed be used to constrain the properties 
of high-z quasars. 

Figure 3d: Based on their recen t study of local 
(U)LIRGs, Gracia-Carpio et al. ( 2006) suggest that the 
^hcn/^hco+ ratio correlates with -Lfir; their results in- 
dicate that HCN may not be an unbiased tracer of star 
formation. We show this relation including their sam- 
ple together with the spiral/star burst sample described 
above and the Cloverleaf. We find no evidence for any 
correlation between £hcn/^hco+ anc ^ ^fir over the in- 
creased luminosity and redshift range. We note that i f 
the relation suggested by Gracia-Carpio et al. (2006) 
were to hold, we would not have been able to detect the 
HCO + (J=1^0) line in the Cloverleaf, as iHCN/-^HCO+ 
would be ~3. 

As discussed above, the £hcn/^hco+ ra ti° is consis- 
tent with unity over a large range in far-IR luminosities. 
Together with multi-transition studies of both molecules 
availa ble for some of the local galaxies (e.g. Seaquist & 
Fraver l2"000|) . this result suggests that HCO + does not 
require special conditions to be excited. Also, it indi- 
cates that HCO + and HCN trace physically similar re- 
gions. Although being consistent with unity within the 
error bars, the ratio of HCN and HCO + luminosities may 
be larger than 1 in the Cloverleaf, and several mecha- 
nisms to explain such a difference have been discussed in 



the literature. It has been argued that HCN emission in 
molecular clouds can be enhanced by mid-IR pumping of 
a 14 ^m vibrational band (Aalto et al. ll995|) . but HCO+ 
can be mid-IR-pumped under very similar conditions 
via a 12 /im vibrationa l band (Gracia-Carpio et al.|2006). 
Also, Gao & Solomon (.20043) have found that this mech- 
anism does not appear to play a major role (their Fig. 5). 
It has also been suggested that the chemical abundance 
of HCN can be enhanced relative to CO and HCO + in the 
ambient X-ray radiation field of a strong AGN (Lepp & 
Dalgarno 1996). It is also possible that the HCO + abun- 
dance is decreased due to the ionizin g field produced by 
cosmic rays (Seaquist & FrayerHoOQ) : Cosmic rays ionize 
H2, leading to the production of H^, which reacts with 
CO to form HCO+. The abundance of HCO+ is thus af- 
fected by ratio of cosmic-ray ionization rate and gas den- 
sity. While a higher ionizing flux favours the production 
of HCO + , it also increases the number of free electrons, 
which leads to a higher probability for recombination: At 
a gas density of 3 x 10 4 cm -3 , an ionizing field compara- 
ble to that of the Galaxy may already be strong enough 
to significantly decrease the abundance of HCO + due to 
disso ciative recombination of H3" (e.g. Phillips & Lazio 
1995), and the ionizing field in the Cloverleaf is probably 
much stronger than the Galactic one. 

Our finding that the HCN-to-HCO + luminosity ratio 
in the Cloverleaf is close to unity implies that the pro- 
cesses discussed here likely do not play a dominant role. 
All these chemical arguments, assume that the HCN and 
HCO + J = 1 — > opacities are low and that the ob- 
served line intensities scale with the underlying molecular 
abundances. Observations of the 13 C baring isotophome- 
res of HCO + and HCN in nearby starburst galaxies have 
shown that the HCO+/H 13 CO+ and HCN/H 13 CN ratios 
are si milar to those seen in C O/ 13 CO (Nguyen-Q-Rieu 
et al. Il992t Wang et al. 120041) . From this it has been 
concluded that that the HCO + and HCN opacities are 
similar to those in CO with t ~ 3 — 4 (see also Henkel 
et al. 119931) . i.e. the emission is optically thick. This 
would imply that ^hcn/^hco+ * s s °l er y determined by 
the relative area filling factors and excitation temper- 
atures of both molecules. As the HCO + (J=1^0) and 
HCN(J = 1^0) excitation is likely to be close to ther- 
malized for densities of n{B.2) ^ 10 5 cm -3 , this would 
naturally explain £hcn/-^hco+ ~ •"■ assuming that both 
molecules trace regions of similar density (i.e., similar 
volume). We conclude that HCO + compares favourably 
with HCN in terms of being a good tracer for dense 
molecular gas even in FIR-bright objects at high red- 
shift. 
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RIGHT ASCENSION (J2000) 

Fig. 1. — VLA detection of HCO"*~(J=l— »0) towards the Cloverleaf quasar at a resolution of 4.0 "x 3.0" (as indicated in the bottom 
left corner). Continuum emission was subtracted. The source is marginally resolved. The cross indicates the geometrical center of the 
CO emission in the Cloverleaf (Alloin et al. 119971 see text). This continuum-subtracted map is integrated over the central 411kms — 1 
(34.375 MHz). Contours are shown at (-3, -2, 2, 3, 4, 5, 6, 7, 8)xct (la = 16/iJy beam -1 ). 




14 15 47.0 46.0 

Fig. 2. — Channel maps of the HCO + (J=l— >0) emission (same region is shown as in Fig. 1). One channel width is 6.25 MHz, or 75kms — 1 
(frequencies increase with channel number and are shown at 25057.228, 25063.478, 25069.728 and 25075.978 MHz). Contours are shown at 
(-3, -2, 2, 3, 4)x<t (lo- = 50/iJybeam- 1 ). The beam size (4.0"x3.0") is shown in the bottom left corner; the cross indicates the same 
position as in Fig. 1. 
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Fig. 3. — HC O+ lum inosity relations for a sample of low-z spi ral an d starburst galaxies (Nguyen-Q-Rieu et al. ll992t Imanishi et al. 120041 : 
Gao fc Solomon 2004a), low— z (U)LIRGs (Gracia-Carpio et al. I2006t Imanishi et al. 20091 . and the Cloverleaf (this work; Solomon et al. 
2003; Weiss et al. 200if). The Cloverleaf luminosities are corrected for gravitational lensing (//l = 11). The solid lines are least squares fits 
to all data except the Cloverleaf. The error bars shown for the Cloverleaf indicate the statistical errors of the line luminosity measurements. 
See text for more details. 



TABLE 1 

CO, HCN, AND HCO+ LINE LUMINOSITIES IN THE CLOVERLEAF. 





Su 

[^Jy] 


L' 

[10 9 K kms-ipc 2 ] 


Ref. 


HCO+(J=1^0) 


193 ± 28 


3.5 ± 0.3 


1 


HCN(J=1^0) 


240 ± 40 


4.3 ± 0.5 


2 


12 CO(J=3^2) 


30000 ± 1700 


42 ± 7 


3 



References. — [1] This work, [2] Solomon ctal. 12003ft . [3] Weiss et al. 120031) . 

Note. — Luminosities arc derived as described by Solomon etal. (1992b): h' [K km s _1 pc 2 ] — 3.25 x 10 7 X I X x x (1 + z)~ 3 , where I 

is the velocity— integrated line flux in Jy kms -1 , Dj_, is the luminosity distance in Mpc (z — 2.55784, Weiss et al. I2003D . and ^ bs ls the observed 
frequency in GHz. All given luminosities arc corrected for this lensing maginification. The HCN and HCO^ luminosities arc corrected for the finite 
source size relative to the synthesized VLA beam (see text). 



